Propagation of nuclear data uncertainties for ELECTRA burn-up calculations 
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The European Lead-Cooled Training Reactor (ELECTRA) has been proposed as a training reactor 
for fast systems within the Swedish nuclear program. It is a low-power fast reactor cooled by pure 
liquid lead. In this work, we propagate the uncertainties in ^^^Pu transport data to uncertainties 
in the fuel inventory of ELECTRA during the reactor life using the Total Monte Carlo approach 
(TMC). Within the TENDL project the nuclear models input parameters were randomized within 
their uncertainties and 740 ^^'^ Pu nuclear data libraries were generated. These libraries are used as 
inputs to reactor codes, in our case SERPENT, to perform uncertainty analysis of nuclear reactor 
inventory during burn-up. The uncertainty in the inventory determines uncertainties in: the long- 
term radio-toxicity, the decay heat, the evolution of reactivity parameters, gas pressure and volatile 
fission product content. In this work, a methodology called fast TMC is utilized, which reduces the 
overall calculation time. The uncertainty in the long-term radiotoxicity, decay heat, gas pressure 
and volatile fission products were found to be insignificant. However, the uncertainty of some minor 
actinides were observed to be rather large and therefore their impact on multiple recycling should be 
investigated further. It was also found that, criticality benchmarks can be used to reduce inventory 
uncertainties due to nuclear data. Further studies are needed to include fission yield uncertainties, 
more isotopes, and a larger set of benchmarks. 
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I. INTRODUCTION 

Accurate and reliable nuclear data (ND) and their un- 
certainties are important in the design, modeling and de- 
velopment of GEN-IV reactors. Out of the six nuclear 
reactors selected by the GEN-IV International Forum, 
three are fast spectrum reactors which have an added 
advantage for transmutation of minor actinides such as, 
americium and curium[l|. The transmutation perfor- 
mances are significantly affected by nuclear data and it is 
important to quantify the resulting uncertainties on the 
integral parameters on fast reactor systems since their 
fuel can be heavily loaded vifith minor actinides y. Dur- 
ing reactor operation, many nuclides are produced as a 
result of different types of neutron-nuclear interactions. 
The long-term changes in the properties of a nuclear re- 
actor core over its lifetime can be determined by these 
changes in fuel composition due to fuel burnup. Knowl- 
edge of isotopic composition during reactor life time is 
important for the estimation of long-term radiotoxicity 
and decay heat of discharged fuel, the evolution of reac- 
tivity parameters, build up of gas pressure, and inventory 
of volatile fission product content (which constitute the 



source term in accident scenarios). For transport and 
storage of spent fuel both an accurate prediction of the 
total decay heat is essential, as well as the uncertainty in 
keff. These are the motivations to investigate the impact 
of nuclear data uncertainties in burnup for the European 
Lead Cooled Training reactor (ELECTRA) [^ , which is a 
low power (0.5 MW) reactor intended for education and 
training purposes. To allow for a small core with cooling 
accomplished by natural circulation an inert matrix zirco- 
nium, plutonium nitride fuel has been chosen for ELEC- 
TRA. The fuel is composed of 60 mol% ZrN and 40 mol% 
of PuN with a TRU vector (238-242 p^^ 24i^^. 35/ g^^g 

/ 23.8 / 11.7 / 7.9 / 1.2 [at. %]), which resembles typical 
spent fuel available at the Swedish interim storage [^. 
The Total Monte Carlo method Q was applied to study 
the effects of uncertainties in "^^^Pu transport and acti- 
vation data on the ELECTRA inventory during burn-up. 
In order to study the effects of transport nuclear data 
uncertainties 740 ^^Spu files from the TENDL project [i| 
have been used. These files include information on trans- 
port data, i.e. data necessary for the calculating neutron 
transport in steady state (cross-sections, angular distri- 
butions, differential data, emission spectra and nu-bar). 
The same cross-section data was used as activation data 
in the depletion calculations. 
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II. TOTAL MONTE CARLO APPROACH 

In the Total Monte Carlo approach [j| uncertainties 
in nuclear model parameter are derived from theoreti- 
cal considerations. They are randomly sampled multiple 
times and each set of parameters are used in a suit of nu- 
clear data codes: TALYS, TANES, TARES, TAPIS Q. 
The output is a large set of random nuclear data files 
which are processed into ENDF format using the TEFAL 
code|j]. These can then be compared with experimental 
data and some are discarded if the deviation from ex- 
perimental data is too large. The accepted random files 
are processed using the ACER module in NJOY, which 
converts the random ENDF nuclear data files into ACE 
files. Theses random files can then be used in neutron 
transport and activation codes such as SERPENT [0] or 
MCNP. The variation of the nuclear data causes distribu- 
tions in the obtained quantities such as fee//, inventory, 
temperature feedback coefficients, and kinetic parameters 
etc. [1|. The spread in these distributions {aobs) is both 
due to variations in ND in the different random files and 
due to the statistical uncertainties from the Monte Carlo 
transport calculation: 



<^obs = <^ND + f^.s 



_ _ tics (1) 

Since the transport code calculates astat, it is possible 
to infer (Tnd, which is an estimated of the uncertainty 
due to ND in the investigated quantity. 



III. SIMULATIONS 

Burnup calculations were performed using the inte- 
grated neutron transport and depletion code SERPENT 
1.1.17 [i]. One advantage of using SERPENT is that no 
distinction is made between the transport nuclear data 
used for calculating e.g. the flux spectrum and the acti- 
vation nuclear data (cross-sections) used for solving the 
Bateman equation. Therefore the same random files are 
used both in the transport part and depletion part of 
the code. Consequently, full information on the corre- 
lations in nuclear data uncertainty propagation between 
these two steps is obtained. As mentioned earlier 740 
^^^Pu files were obtained from the TENDL project and 
processed into ACE files with Njoy99.336 code. JEFF3.1 
general purpose library was used for all the other iso- 
topes. The other input data used in this work are the 
decay and fission yield data read from the JEFF-3.1 nu- 
clear data library. These files were not varied in this work. 
The calculations were performed for 20 inactive and 50 
active cycles with 5000 neutrons/cycle. Burnup calcula- 
tions were performed by invoking the predictor-corrector 
module in SERPENT for eight burnup steps from 0-55 
Mwd/kg corresponding to the beginning of life (BOL) 
and end of life (EOL) of the reactor respectively Q. For 
decay calculations, eighteen steps were used from 0—10^ 
years after the reactor EOL. Continuous full power op- 
eration is assumed, up till the end of life (EOL) of the 



reactor. The EOL is reached when the reactivity drums 
has reached their end position. This was calculated to 
be at 55 at MWd/kg in ref. 3, which occurs after 7442 
full power days. In the simulations performed for this pa- 
per the drums were not rotated but the simulation were 
executed with the drums in their start position. 



IV. FAST TMC IN BURNUP CALCULATIONS 

In ref. Q fast TMC methodologies are described in 
detail and in this section only the burnup aspect will be 
discussed. From Eq. (1) we see that we need an estimate 
on (7 Stat, but statistical uncertainties are not propagated 
in burnup codes. Therefore we run our transport code in 
two sets with the only difference that in the first set the 
code is run Nl times with constant ND and in the second 
set the code is run the N2 times, each time with a different 
random file (N2 = the number of random files). In our 
case Nl= 700 and N2 = 740. FIG. 1 shows the result, 
where the set with constant nuclear data has a narrower 
distribution. Its spread is only due to statistics in the 
calculation, <Tstat- The set that correspond to the 740 
random files has a wider distribution and its spread, aots, 
is both due to statistics and ND as outlined in Eq.(l). 
With this information ajyo can be inferred. Since we 
have many random files, the uncertainties on the spreads 
{Aastat and Across) are small. As a consequence, we can 
reduce the statistical requirement on each run, lowering 
the number of neutrons, hence reducing the computation 
time. This is why the methodology is referred to as fast 
TMC. In this work the full core burnup calculation took 
24 h on 60 cores. 



■ 740 random Pu-23<i 

740 random Pu-23? 

□Constat ND 
— Constant ND 




FIG. 1. The distribution in fce// in the two sets of simula- 
tions. The blue bars represent the 740 random files, where 
each simulation has unique ^^'^ Pu ND. The light bars repre- 
sent the outcome from the 700 simulations where the ND was 
kept constant. All simulations have unique seeds. 



It is also of interest to know Actjvd, i.e., the uncer- 
tainty on the uncertainty due to nuclear data. We first 
compute the uncertainty of the variance due to nuclear 
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data (V/vd): 



(AVnd) = [AVobsf + iAV,i 



(2) 



Where Vobs is the variance in the random files, Vstat is 
the variance for the set with constant nuclear data and A 
is the associated uncertainties. The uncertainty in Vobs 

is: 



S„ 



AVobs = V2^ 



(3) 



Where capital S is the sample variance. The uncer- 
tainty of Vstat is: 



AK. 



V2. 



^stat 



Combining the above equations: 

2 




V2 



^stat 



(4) 



(5) 



Normally, what is asked for is the uncertainty in uncer- 
tainty due to ND , which is: 



Acr^£, 



AV, 



ND 
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N2 



Ni 



2cr 



ND 
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■^obs' 



'ND 



(6) 



As can be observed from the equation it is more im- 
portant to have many random files than to improve the 
statistical spread in each run. The number of random 
files should normally be higher than the number of runs 
with constant ND. 



V. RESULTS AND DISCUSSION 

In previous work on ELECTRA, [^ it was observed 
that the total reactivity loss at EOL was 9800 ± 200 pcm 
with 13.8% of the initial plutonium converted into fission 
products or americium. In FIG. 2, we present the main 
results on the uncertainty in inventory and uncertainty in 
decay heat. When we look at the evolution during bur- 
nup it is seen that there is only small uncertainties for 
major actinides and consequently there is only a small 
uncertainty in the radiotoxicity. For the minor actinides 
larger uncertainties are observed e.g. 11% in ^"'^C/. This 
has little effect on integral parameters such as radiotoxic- 
ity (not shown) or decay heat (FIG. 2, top). However the 
large uncertainties in the minor actinides due to '^^^Pu 
data should be investigated in more detail since it might 
have effects in multiple recycling scenarios, both on de- 
cay heat and radiotoxicity as well as on different safety 
parameters such as the Doppler coefficient. 




rl 




^__x 


\ 






Ss 10-' 




c -^ 




---Actinides 








^ S 




Fission products 








10 




-^plutonium 
— total 
















„ 






y 






m 


^ y^ / 








/ / - 


t: 

a 




/ 


/ / 










S 10 ' 


// 














I 




"•= 2 


y .y 






, 


Jio' 


/■'X 


--Pu- 


239 


-g 


y .y 


— Pu- 


240 






/■/ 


-^Am 


-241 




S'o' 


yy/ 


^Am 


-243 
-244 


1 


5 


■y 







Time from reactor start [YEARS] 

FIG. 2. (Top) relative decay heat uncertainties due to trans- 
port (and activation) ^^^Pu data (%) as a function of time 
(years). (Bottom) relative uncertainties in the inventory due 
to transport (and activation) ^^^^ Pu data for some selected 
isotopes as a function of time (years). At 20 years the reactor 
reaches its EOL. 



A. Benchmarks for inventory 

In ref. Q benchmarks are used to select the best ran- 
dom file and in ref. [8|, a new method to reduce the 
uncertainty in k^ff due to nuclear data using criticality 
benchmarks is presented. In this work we use the same 
methodology to investigate if criticality benchmarks can 
be used to reduce uncertainties in the inventories for the 
different isotopes. The Jezebel benchmark Q was used as 
the criticality benchmark. All random files were also run 
on that benchmark. To use the benchmark to reduce the 
uncertainty, a correlation between the benchmark and the 
investigated parameter (in this case the number density) 
has to be found. FIG. 3 shows examples of correlations 
between the number density in ^'^^Am and '^^^Pu and 
the keff in Jezebel for the different benchmarks. A cor- 
relation exists between the ^'^^Am number density and 
Jezebel fee//- Consequently, the Jezebel Benchmark can 
be used to reduce uncertainties in the '^'^^Am inventory. 
This is not the case for the '^^^Pu inventory. 

For each random file a goodness number is calculated 



X 



2. 






(7) 



Where Cj is the calculated fee// for the the random 
file J, and E is the nominal value for the benchmark, in 
this case 1.0000. An acceptance criteria based on x^ can 
be used, which is illustrated in FIG. 4, where it is shown 
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FIG. 3. The ^^'^Am (top) and '^^^Pu (bottom) inventories for 
the 740 ^^^Pu random files are plotted (each blue circle rep- 
resenting one random file) against the calculated k^ff in the 
Jezebel Benchmark for the same files. The Jezebel benchmark 
has a nominal value of 1.000. 
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FIG. 4. The uncertainty in "^"^^ Am inventory due to nuclear 
data vs. the maximum allowed x for the random file on the 
Jezebel benchmark. 



that with a more rigid acceptance criteria the uncertainty 
due to ND is reduced. It is found that a 25% reduction 
in uncertainty is achieved for "^^^Am, however no signifi- 
cant reduction in uncertainty was achieved for ^^^Pw. All 
actinides were tested for correlations in all burnup steps 
and in most cases, except for "^^^Pu and ^^^Pu, correla- 
tions were observed. It was found that the correlation 
was quite constant with burnup. 

VI. CONCLUSIONS AND OUTLOOK 

Using SERPENT, fast TMC has successfully been per- 
formed for burnup in a full core 3-D model of ELECTRA 
using 740 ^^^Pu transport-data random files. This was 
shown to be achievable with quite limited computer re- 
sources, showing the strength of fast TMC. In the work 
the same random files were used in the depletion calcula- 
tion and the transport calculation. It was found that the 
uncertainty in the long-term radiotoxicity, decay heat, 
gas pressure and volatile fission products is small for 
^■^^Pu transport data uncertainties. The uncertainties on 
some minor actinides are large enough that their impact 
on safety parameters in multiple recycling should be in- 
vestigated further. Decay data and fission yield data was 
not varied. A next step would be to include fission yield 
uncertainties by e.g. combining TALYS and GEF codes. 
Criticality benchmarks can be used to select acceptable 
ND files and thereby reduce inventory uncertainties. In 
a next step a larger set of benchmarks should be used. 
It was found that the Jezebel benchmark could be used 
to reduce the uncertainty in inventory for some isotopes, 
but not for all. The reason for this has to be investigated 
further. 
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